par(mfrow=c(1, 2))
data<-read.csv('data/时间序列分析——基于R（第2版）案例数据/csv/A1_7.csv')

# 设置时间序列变量
number<-ts(data$number, start=1900)
# 绘制时序图
plot(number, type="o", pch=5, col='#39CBB4', main="时序图", xlab="year", ylab="percent")

# 绘制自相关图
acf(number, lag.max=20)$acf

# 纯随机性检验
for( k in 1:3) print(Box.test(number, lag=6*k, type="Ljung-Box"))

